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INTRODUCTION 

The emphasis of this grant was focused on precision ephemerides for the Global Positioning 
System (GPS) satellites for geodynamics applications. During the period of this grant, major 
activities were in the areas of thermal force modeling, numerical integration accuracy improvement 
for eclipsing satellites, analysis of GIG ’91 campaign data and the Southwest Pacific campaign data 
analysis. 

THERMAL FORCE AND ECLIPSING ORBITS 

Papers on thermal imbalance force modeling and eclipsing orbit analysis were presented at the 
AIAA Space Flight Mechanics Conference held in February 1992 at Houston, Texas. Details are 
included in Appendix VII and Appendix VIII of the Final Report for NASA Grant No. NAG5-940. 

GIG ’91 DATA ANALYSIS 

An intensive global GPS observation campaign called "The First GPS IERS and Geodynamics 
Experiment - GIG ’91" was organized by the Jet Propulsion Laboratory for a period of three weeks 
during early 1992. The campaign spanned part of GPS Weeks 576, 577, 578 and part of Week 579 
starting from January 22, 1992, to February 13, 1992. About 120 stations distributed globally, 
collected data during this period using a variety of GPS receivers which included Rogues, TI-4 100’s, 
Trimbles, Ashtechs, Mini Mac 2816 AT’s, WM-102’s; in addition, data from the five continuously 
operating DMA stations were also available. From among these, data from 20 globally distributed 
stations were chosen for orbit determination and Earth orientation parameter determination 
experiments at UT/CSR. Of these, 17 sites had Rogue receivers, two TI-4100 and one Mini Mac 
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2816AT. Orbit determination experiments included computation of orbits using data from a reduced 
set of global stations and from an expanded set, and comparing the bsaelines and other criteria to 
assess the orbit accuracy. In addition to the GIG’91 data set mentioned above, data collected during 
the 1989 South West Pacific Campaign also were used. Results of this study were presented at the 
1991 AGU Spring Meeting held at Baltimore, Maryland, during May 1991. Further details are 
included in the Appendix. 

Using the same GIG’91 data set several experiments were performed at CSR to determine the 
pole positions ( X p , Y p ) during this time period. The main objective of the experiments was to 
explore strategies for determining the Earth orientation parameters (EOP), polar motion and UT1, 
using global GPS data. Short arc and long arc with sub arc parameters were considered as possible 
approaches. The estimated pole positions were compared with SLR and VLBI determined polar 
motion series in order to assess the quality of these determinations. Results of this study were 
presented in the special workshop held in Ahrweiler, Germany, during August 1991 and at the XX 
General Assembly of the IUGG during August 1991 in Vienna, Austria. A copy of the presentation 
is included in the Appendix. 

SWP CAMPAIGN DATA ANALYSIS 

Analysis of data collected during the South West Pacific (SWP) campaigns contributed to 
accuracy assessment of the GPS satellite orbits. The SWP campaign data were collected during the 
summers of 1988, 1989 and 1990 and these data have been analyzed in various stages at CSR. These 
data sets complement the CIGNET data in forming a better global distribution and facilitates various 
orbit and baseline experiments. Studying the repeatability of baselines between various sites in the 
SWP network provides one method to evaluate the accuracy level of the computed orbits. In depth 
analysis of few weeks of data collected at the sites in Tonga region (SWP) during the above three 
years were performed at CSR and the results indicate that the computed GPS orbits yield baseline 
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repeatability at the 10 to 20 parts per billion (ppb) level. A brief summary of these results were 
presented at the 1991 AGU Fall Meeting held in San Francisco during December 1991. Copy of the 
presentation is included in Appendix. 

CONCLUSION 

Based on the research efforts at CSR in the area of precise ephemerides for GPS satellites, the 
following observations can be made pertaining to the status and future work needed regarding orbit 
accuracy. There are several aspects which need to be addressed in discussing determination of 
precise orbits, such as force models, kinematic models, measurement models, data 
reduction/estimation methods etc. Although each one of these aspects has been studied at CSR in 
research efforts under this (and the previous) grant, only points pertaining to the force modeling 
aspect are addressed here. 

Dynamic Modeling - Current Status 

At present the following known forces are modeled in routine computation of GPS satellite 
orbits: 

(1) Nonspherical Earth gravity acceleration represented by one of the extant gravity fields such 
as GEM-L2, GEM-T1 or GEM-T2 truncated to 8x8 is adequate. However, increasing the degree 
and order and/or tuning the field (at least the resonance coefficients) for GPS orbits may slightly 
improve the accuracy of this perturbation modeling. 

(2) represented as point masses are adequate. 

(3) force modeling and the geometric tide effects must be considered in the measurement 
models. 

(4) Perturbation due to solar radiation pressure is adequately modeled by the ROCK4 models; 
however, it is necessary to scale these accelerations by at least one adjustable parameter. 
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(5) Perturbation due to venting of the heat source in the GPS satellite is modeled as ’Y-bias 
acceleration’ scaled by an adjustable parameter. 

Values of weekly estimates of this scale parameter considered over a period of time do not 
exhibit any systematic trend except for a somewhat weak correlation with the eclipsing period. 
Experience shows that it is necessary to include this perturbation along with an adjustable parameter 
in order to obtain good fit of the data. However, in a long arc solution, considering several sub arcs 
for this perturbation does not significantly improve the rms of fit. 

Following are some of the factors which could contribute to additional improvement of orbit 
accuracy beyond the current level. 

(1) Perturbation due to imbalance in thermal radiation has been shown to cause differences in 
orbit prediction at the level of a few meters over a period of about one week or more. Although 
small, inclusion of this perturbation may help in achieving baseline accuracies (and/or 
repeatabilities) at parts per billion level. But the difficulty in considering this perturbation routinely 
is due to the fact that nonlinear partial differential equations (heat equations) must be solved 
simultaneously with the ordinary differential equations of motion, which causes significant 
complications in algorithm and computation even for a modest approximation of the satellite 
configuration. Hence, careful evaluation of the costs and benefits of including this perturbation is 
needed. 

(2) handling discontinuities in function values (occurring in SRP acceleration at shadow 
crossings), can be overcome by the ad hoc modification of the integrator back difference table. 
However, inclusion of this modification did not seem to improve the rms of fit or the prediction error 
in real data processing, although improvements were obvious in simulation studies. The reason for 
this anomaly is not known at present and will have to be investigated before this feature can 
routinely be included in orbit computation. 
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Dynamic Modeling - Future Study 

There are indications (evidenced by discontinuities in daily/weekly solutions and by prediction 
errors) to the effect that all the perturbations described above do not completely or exactly represent 
all the forces acting on the GPS satellites. There may be other unmodeled forces such as 
unintentional thrusting (due to outgassing, momentum dumping, attitude correction etc.) or due to 
other natural phenomena. One of the ways in which such unknown and unmodeled perturbations 
could be accounted for in orbit computation is to estimate empirical accelerations. Such an approach 
needs detailed analysis in the future. 
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APPENDIX 


1. "GPS Ephemeris Accuracy Improvement From Global Data Set," AGU 1991 
Spring Meeting, May 28-31, 1991, Baltimore, Maryland. 

2. "Analysis of GIG 91," presented at GIG Workshop, Ahrweiler, Germany, 
August 1991. 

3. "Global GPS Orbit Determination," XX General Assembly, IUGG, August 
11-24, 1991, Vienna, Austria. 

4. "Geodetic Analysis of GPS Measurements Near the Tonga Trench: 
1988-1990," AGU 1991 Fall Meeting, December 1991, San Francisco, 
California. 
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INCLUDING THE FIVE DMA STATIONS (TOTAL 19) ; GPS WEEK-499 


GLOBAL TRACKING NETWORK DATA FOR WEEK 499 
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SOLUTION CHARACTERISTICS FOR THE GLOBAL DATA SET 
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TROMSO-ONSALA VECTOR BASELINE ACCURACY 
COMPARISON WITH VLBI SOLUTION 
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RMS OF TRAJECTORY DIFFERENCES 
9 GS vs 20 GS SOLUTIONS 
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(Pred vs Fit) in W500 for PRN-03 
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ECF EPH DIFFERENCES (WEEK 499) 
(UT/CSR DD PHASE SOLN - NSWC EPH) 
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GIG' 91 Workshop 

Agenda 


I 

Monday, 5 August 1991 

Arrival of the participants 

19.00 - 20.00 Open-air reception by good weather otherwise in the restau- 
rant of the hotel 


Tuesday, 6 August 1991 

08.00 - 09.00 Registration 

09.00 - 10.30 Opening Session (Chair: H. Seeger, I FAG) 

Greeting and Remarks from Hosts: H. Seeger (IFAG) 

Greeting and Remarks from IERS-Background and objectives of 
the IERS GPS Campaign: M. Feissel? (IGN) 

Greeting and Remarks from the IERS GPS Technique 
Coordinator: W. Melbourne (JPL) 

IERS Latest UTPM Results for the GIG Experiment Period: 

M. Feissel? (IGN) 

10.30 - 11.00 Coffee break 

11.00 - 12.30 The Field Campaign (Ch: W. Schlilter, IFAG/R. Neilan, JPL) 

Summary of the Campaign-Planning, Standards, Procedures, 
Operations, Performance, etc: R. Neilan/S. Fisher (JPL) 

Reports from Various Operations Teams (Representatives from 
organizations that participated in field ops will be asked to 
give a brief account of their experiences, insights, etc.) 

12.30 - 14.00 Lunch 

14.00 - 15.30 The Field Campaign (cont.) 

Continued Reports from Ops Teams 

Discussion: Lessons Learned, Recommendations for Future 
Campaigns (Standards, Procedures, Documentation, etc.) 

15.30 - 16.00 Coffee break 

16.00 - 18.00 Data Analysis - I: Preparation 

Status of Data Pre-Processing & Distrib., Site Ties 
(S. Fisher) 

Reports by Pre-Processing Centers (All) 

Summary of Analysis Standards Proposed by JPL (G. Blewitt) 
Discussion of above and related topics 

Open-air grillparty 


19.00 



08.30 - 10.30 Data Analysis - II: Results - The plan is to organize this all-day session by ^ 
topics. Those with results in several areas will therefore give several short 
presentations. The main topics are: 

1. Descriptions of software and modeling strategies 

2. Baseline Solutions 

3. GPS Orbit Solutions 

4. Earth Orientation Solutions 

5. Geocenter Solutions 

6. Discussion, Ccnparison & Analysis of Results 

As of 25 July, the known groups planning to present results are: 

U. Texas - B. Schutz 
DGFI - K. Kaniuth 
NSWC - E. Swift 
UNAVCO - C. Rocken 
MIT - R. King 

JPL - Blewitt, Lichten, Lindqwister, Webb, Yunck 
There is still space for more 


10.30 - 11.00 

11.00 - 12.30 

12.30 - 14.00 

14.00 - 15.30 

15.30 - 16.00 

16.00 - 18.00 

19.30 

20.00 - 22.30 

22.30 


Coffee break 

Data Analysis - II: Results (cent.) 

(See Above) 

Lunch 

Data Analysis - II: Results (cont.) 

(See Above) 

Coffee break 

Data Analysis - II: Results (cont.) 

(See Above) 

Bus departure for wine-tasting 
Wine-tasting in Ahrweiler 
Busses return to hotel 


Thursday, 8 August 1991 


08.30 - 10.00. Special Topics (Preliminary list): 


Multipath studies with GIG data (C. Rocken) 

Do high latitude sites pose special problems? (TBD) 
triplications of SA/AS for global GPS measurements (T. Yunck) 
Status of 3CAR GPS geodesy project in Antarctica (J. Manning) 
Sirmary of early GIG results (TBD) 

10.00 - 10.30 Coffee break 

10.30 - 12.30 Open Discussion (Suggested topics): 


Implications of early GIG results 

Establishing a GPS-based global reference frame Integrating GPS into IERS 
v - }•>•. operations and products Plans for a follow-up GIG" 91 workshop before IGS 

campaign? 

12.30 ' Lunch 


14.30 


Sightseeing tour through the Ahr-Valley 


ANALYSIS OF GIG 91 
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Pole position error estimate (RMS): 1 mas 



GIG 91 DATA 
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GLOBAL GPS ORBIT DETERMINATION 
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SOFTWARE VALIDATION 
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SOFTWARE VALIDATION (continued) 
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GEODETIC PROCESSING (COMPUTER: CRAY Y-MP) 
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DATA PREPROCESSING (COMPUTER: VAX) 
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GIG '91 

GPS Experiment for IERS and Geodynamics 
ROGUE RECEIVER LOCATIONS 



ADDITIONAL ROGUES IN CALIFORNIA PERMANENT ARRAY: 
PINON, SCRIPPS, JPL MESA 




GIG 91 DATA 
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GPS ANALYSIS MODELS 



Empirical Forces: 

• Along track (CT) 

• Radial (once/rev) 

• Normal (once/rev) 



GPS MODELS (continued) 
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REPEATABILITY BASED ON 4 DAYS 

8 STA, 28 BASELINES x vert(M) 
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Polar Motion Comparison During the GIG Campaign 











EODETIC ANALYSIS OF GPS 
IEASUREMENTS NEAR THE 
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DECEMBER 1991 
SAN FRANCISCO 




| continental crust A A A subduction zone spreading ridge offset by transform faults 

• •volcanoes uncertain plate boundary - — — collision zone 

v earthquake zone ► movement of plate 21 Plate boundaries and active zones of the Earth's crust 












SOUTHWEST PACIFIC DATA CHARACTERISTICS 
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SWP Network: Trimble SST (11 sites) 

North American Evaluation Network: Trimble SST & TI-4100 (3 sites) 
SA on, 8 day campaign, 13 satellites 



SOUTHWEST PACIFIC PROCESSING STATUS 
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Figure 2. 1988 Southwest Pacific Campaign Global Sites 


SWP-90 ANALYSIS STRATEGY AND EXPERIMENTS 





SWP-90 SOLUTIONS 
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Precision estimate: 

• From above 3 lines: 10-20 parts per billion 

• From 103 lines in the global network: 2.3 cm + 8 ppb*length 



SWP-90 EVALUATION SOLUTIONS 



Accuracy estimate: 10-30 parts per billion 

• Receiver mix used in evaluation 

• Evaluation sites include mix of SLR and VLB I 



SWP-88 SOLUTIONS 
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Repeatabilities (combined 444 and 445): 

• Rarotonga/Tongatapu: 4.6 cm 

• Rarotonga/Vava’u: 4.1 cm 

• Rarotonga/W. Samoa: 6.9 cm 
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